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ABSTRACT 

The Extragalactic Background Light (EBL) includes photons with wave- 
lengths from ultraviolet to infrared, which are effective at attenuating gamma 
rays with energy above ~ 10 GeV during propagation from sources at cosmo log- 
ical distances. This results in a redshift- and energy-dependent attenuation of 
the 7-ray flux of extragalactic sources such as blazars and Gamma-Ray Bursts 
(GRBs). The Large Area Telescope onboard Fermi detects a sample of 7-ray 
blazars with redshift up to 2; ~ 3, and GRBs with redshift up to 2 ~ 4.3. Using 
photons above 10 GeV collected by Fermi over more than one year of observations 
for these sources, we investigate the effect of 7-ray flux attenuation by the EBL. 
We place upper limits on the 7-ray opacity of the Universe at various energies 
and redshifts, and compare this with predictions from well-known EBL models. 
We find that an EBL intensity in the optical-ultraviolet w avelengths as great as 



predicted by the "baseline" model of IStecker et al.l ( l2006l ) can be ruled out with 
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high confidence. 

Subject headings: diffuse radiation - dust, extinction - Gamma rays: general 
Opacity - Galaxies: active - Gamma-ray burst: general 



1. Introduction 

The Fermi Gamma Ray Space Telescope was launched 2008 June 11, to provide an un- 
precedented view of the 7-ray Universe. The main instrument onboard F ermi, the Large Area 



Telescope (LAT), offers a broader bandpass (~ 20 MeV to over 300 GeV: lAtwood et al.ll2009l ) 



and improved sensitivity (by greater than an order of magnitude) tha n that of its predecesso r 



instrument EGRET onboard the Compton Gam ma Ray Observatory (JThompson et al.lll993l ) 



and the Italian Space Agency satellite AGILE (iTavani et al.ll2008l ). which was launched in 
2007. The LAT observes the full sky every 3 hr in survey mode leading to a broadly uniform 
exposure with less than ~ 15% variation. The Gamma-ray Burst Monitor, the lower energy 
(~ 8 keV - 40 MeV) instrument onboard Fermi, observes the full un-occulted sky at all 
times and provides alerts for transient sources such as GRBs. 

A major science goal of Fermi is to probe the opacity of the Universe to high-energy (HE) 
7-rays as they propagate from their sources to Earth. Such energetic photons are subject to 
absorption by production of e lectron-positron (e~e'^) pairs while inte r acting with low energy 



cosmic background photons (JNishikovl Il96ll : iGould fc Shrederl Il966l : iFazio fc Steckerl Il970l ) 



if above the interaction threshold: ethr = {2meC^)'^ / {2E{1 — /i)) where e and E denote the 
energies of the background photon and 7 ray, respectively, in the comoving frame of the 
interaction, rrieC^ is the rest mass electron energy, and 6 = arccos(/i) the interaction angle. 
Because of the sharply peaked cross section close to threshold, most interactions are centered 
around e* ~ 0.8(i?/TeV)~^eV for a smooth broadband spectrum. Thus, the extragalactic 
background light (EBL) at UV through optical wavelengths constitutes the main source of 
opacity for 7-rays from extragalactic sources (Active Galactic Nuclei: AGN and GRBs) in 
the LAT energy range. The effect of absorption of HE 7-rays is then reflected in an energy- 
and redshift dependent softening of the observed spectrum from a distant 7-ray source. The 
observation, or absence, of such spectral features at HEs, for a source at redshift z can be 
used to constrain the 77 — )■ e~^e~ pair production optical depth, T^-y{E, z). 

The EBL is dominated by radiation from stars, directly from their surface and via 
reprocessing by dust in their host galaxies, that accumulated over cosmological evolution. 
Knowledge of its intensity variation with time would probe models of galaxy and star forma- 
tion. The intensity of the EBL from the near-IR to ultraviolet is thought to be dominated by 



-7- 



direct starlight emission out to large redsliifts, and to a lesser extent by optically bright AGN. 
At longer wavelengths the infrared background is produced by thermal radiation from dust 
whi ch is heated by sta rlight, and also emission from polycyclic aromatic hydrocarbons (see 



e.g. 



Driver et al.l 120081 ) 



Direct measureme nts of the EBL is difficu lt due to contamination by foreground zodiacal 
and Galactic light (e.g.. lHauser fc Dwekll200ll ). and galaxy counts result in a l ower limit since 
the number of unresolved sources is unknown (e.g.. lMadau fc Pozzettil 120001 ). Furthermore, 
evolution of the EBL density in the past epochs {z > 0) that is required to calculate the 
7-ray flux attenuation from distant sources cannot be addressed by measuring the EBL 
density at the present epoch (z=0). Hence, several approaches have been developed to 
calculate the EBL density as a function of redshift. The models encompass different degrees 
of complexity, observational constraints and data inputs. Unfortunately, the available direct 
EBL measurements do not constrain these models strongly at optical-UV wavelengths due 
to the large scatter in the data points. A description of the different models is beyond the 
scope of this work; we refer the reader to the original works on th e various EBL models (e.g. 



Salamon fc Stecker' 19981: Stecker et all 2006: Kneiske et al. 2002. 2004; Primack et al. 



Gilmore et al.i i2009i : iFranceschini et al.l l2008l : iRazzaque et al.l 



2009 



Finkeet al.l 120101 ). We 



2005 



note that all recent EBL models, and in particular all models used in this paper, use almost 
identical parameters of a ACDM cosmology model. 

For the analyses presented in this work we have made use of the optical depth values 
t{E, z) provided by the authors of these EBL models. These mod els are available via web- 
page s a analytical appr o ximat ions (as in, e.g.. lStecker et al.l (120061 )). published tables (as in. 



e^ 



Franceschini et al.l (l2008l)) or via pri v ate communications (which i s the case fo r 



_e^ 



Salamon fc Stecked fll998h : IPrimack et al.l ( 120051 ): iGilmore et al.l (l2009f l: iFinke et al.l (I2OIOI ) 
for this work). Since the optical depth values are usually available in tabular form, for ex- 
act values of observed energy E and redshift z a linear interpolation of t{E, z) is used for 
arbitrary values of E and z in our calculations below. 

The range of predictions by these EBL models is illustrated in Figure [1] as a function of 
observed 7-ray energy for sources at different redshifts. The Universe is optically thin (r^^ < 
1) to 7-r ays with energy b elow c^ 10 GeV up to redshift 2; ~ 3, independently of the model 
(see also iHartmannl (120071 )). This is due to the rapid extinction of EBL photons shortwards 
of the Lyman limit. Gamma rays below ~ 10 GeV are not attenuated substantially because 
of faint far-UV and X-ray diffuse backgrounds. 



http: //www. physics . adelaide .edu. au/~tkneiske/Tau_data.html for Kneiske 



2004; 



http : //www . phy . ohiou . edu/'f inke/EBL/index . html, for Finke et al. 2010 
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Fig. 1. — Attenuation as a fu nction of observed g amma-ray energy for the EBL models of 
(jFranceschini et al.ll2008[ l and ( iStecker et al.ll2006l ). These models predict the minimum and 
maximum absorption of all models in the literature, and thus illustrate the range of optical 
depths predicted in the Fermi-LAT energy range. 
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The primary sources of HE extragalactic 7-rays are blazars and gamma ray bursts 
(GRBs). Blazars are active galactic nuclei (AGN) with relativistic plasma outflows (jets) 
directed along our line of sight. GRBs are associated with the core collapse of massive stars, 
or might be caused by binary mergers of neutron stars or neutron star - black hole systems. 
Some GRBs produce beamed high-energy radiation similar to the case of blazars but lasting 
for a short period of time. GRBs have not been used to constrain EBL absorption during 
the pie-Fermi era mainly because of a lack of sensitivity to transient objects above 10 GeV. 
The sensitivity of EGRET decreased significantly above 10 GeV, and the field-of-view (FoV) 
of TeV instruments is small (typically 2 — 4°) to catch the prompt phase where most of the 
HE emission occurs. The new energy window (10 — 300 GeV) accessible by Fermi, and the 
wide FoV of the LAT, makes GRBs interesting targets to constrain EBL absorption in this 
energy band. 

Evaluating the ratio of the putatively absorbed to unabsorbed fiuxes from a large number 
of distant blaz ars and GRBs obse rved by Fermi could result in interesting E BL constraints, 



as pr oposed by lChen et al.l ( 120041 ). although intrinsic spe ctral curvatur e (e.g.. lMassaro et al. 



20061 ) or redshift dependent sou rce internal absorption (IReimeii 120071 ) could make this, or 
similar techniques, less effective. iGeorganopoulos et al.l ( l2008l ) have proposed that Compton 
scattering of the EBL by the radio lobes of nearby radio galaxies such as Fornax A could be 
detectable by the Fermi-LAT. If identified as unambiguously originating from such process, 
a LAT detection of Fornax A could constrain the local EBL intensity. 

Because the e-folding cutoff energy, E{t^^ = 1), from 77 pair production in 7-ray 
source spectra decreases with redshift, modern Cherenkov 7-ray telescopes are limited to 
probing EBL absorption at low redshift due to their detec tion energy thresholds typically at 
or below 50 GeV to 100 GeV (JHinton fc HofmannI 120091 ) . Ground-based 7-ray telescopes 
have detected 35 extragalactic sources to datcl, mostly of the high-synchrotron peaked 
(HSP) BL Lacertae objects type. The most distant sources seen from the ground with 
a confirmed redshi ft are the flat spectrum radio qua sar (FSRQ) 3C 27 9 at z = 0.536 
flAlbert et all boosi ) and PKS 1510-089 at z = 0.36 JWagner et al.lboioh . Observations 
of the closest sources at multi-TeV energies have been effective in placing limits on the 
local EBL at mid-IR wavelengths, while spectra of more distant sources generally do not 
extend above 1 TeV , and therefore probe the optical and near-IR starlight peak of the in- 



tervening EBL (e.g..lStecker fc de Jageiill993l: IStanev fc Franceschinilll998l: ISchroedter 



2005 



2007 



Albert et al. 



Aharonian et al. 1999, 2002; Costamante et al. 



2008 



Krennrich et al. 



2008; 



2004J: lAharonian et al.ll2006l : iMazin fc Raue 
Finke fc Razzaquell2009h ~ 



e.g. http : //www . mpi-hd . mpg . de/hf m/HESS/pages/home/sources/ http : //www . mppmu . mpg . de/~rwagner/sour ces/ 
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The starting point for constraining the EBL intensity from observations of TeV 7-rays 
from distant blazars with atmospheric Cherenkov telescopes is the assumption of a reasonable 
intrinsic blazar spectrum, which, in the case of a power law, dN/dE oc E^^'"* for example, 
that is not harder than a pre-specified minimum value, e.g.. Tint > Tmrn = 0.67 or 1.5. 
Upper limits on the EBL intensity are obtained when the reconstructed intrinsic spectral 
index from the observed spectrum, Tobs, presumably softened by EBL absorption of very high 
energy (VHE) 7-rays, is required to not fall below Tint- The minimum value o f F has been 
a matter of much debate, being reasoned to be Tint = 1-5 by lAharonian et al.l (120061 ) from 



simple shock acceleration theory and from the observed spectral energy distribution (SED) 
properties of blazars, while IStecker et al.l (120071 ) argued for harder values (^l ess than 1.5) under 



specifi c conditions based on more detailed shock acceleration simulations. iKatarzyhski et al. 
( I2OO6I ) suggested that a spectral index as hard as Tint = 0.67 was possible in a single-zone 



leptonic model if the underlying el ectron spectrurn respo nsible for inverse- Compton emission 
had a sharp lower-energy cutoff. iBottcher et al.l ( l2008l ) noted that Compton scattering of 
the cosmic microwave b ackground radiation by extended jets could lead to harder observed 
VHE 7-ray spectra, and lAharonian et al.l (120081 ) have argued that internal absorption could, 
in some cases, lead to harder spectra in the TeV range as well. 

A less model dependent approach uses the (unabsorbed) photon index as measured in 
the sub-GeV range as the intrinsic sp ectral slope at GeV -TeV energies. This method has 



recentlv been apphed to PG 1553+113 ( JAbdo et al.ll2010al l and lES 1424+240 ( JAcciari et al. 



2OIOI : iPrandini et al.ll2010[ ) to derive upper limits on their uncertain redshifts, and to search 
for EBL-induced sp ectral softening in Fermi observations of a sample of TeV-selected AGN 
(lAbdo et allboiOil ). 



Attenuation in the spectra of higher redshift objects {z > 1) may be detectable at the 
lower energies that are accessible to the Fermi-LAT, i.e., at -E ~ 10 — 300 GeV. Gamma rays 
at these energies are attenuated mainly by the evolving UV background, which is produced 
primarily by young stellar populations and closely traces the global star-formation rate. 
Observations with Fermi of sources out to high redshift could therefore reveal information 
about the star-formation history of the Universe, as well as the uncertain attenuation of UV 
starlight by dust. 

In this paper we present constraints on the EBL intensity of the Universe derived from 
Fermi-LPCT observations of blazars and GRBs. The highest-energy 7-rays from high redshift 
sources are the most effective probe of the EBL intensity, and consequently a powerful tool 
for investigating possible signatures of EBL absorption. In contrast to ground-based 7-ray 
detectors, Fermi offers the possibility of probing the EBL at high redshifts by the detection 
of AGN at > 10 GeV energies out to z > 3, and additionally by the detection of GRB 
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080916C at a redshift of ~ 4.35 (lAbdo et al. 



2009 



2009al : ICreiner et all 120091 ). GRBs are known 



GRB 090423 is the current record holder 



to exist at even higher redshifts (iTanvir et al. 
with z ~ 8.2). Therefore observations of these sources with Fermi are promising candidates 
for probing the optical-UV EBL at high redshifts that are not currently accessible to ground- 
based (Cherenkov) telescopes. 

In Section |2] we describe our data selections, the Fermi LAT AGN and GRB observations 
during the first year of operation and analysis, and we discuss potential biases in the selection. 
Our methodology and results are presented in Section [31 We discuss implications of our 
results in Section HI and conclude in Section |5l 

In the following, energies are in the observer frame except where noted otherwise. 



Observations and data selection 



The Fermi LAT is a pair-conversion detector sensitive to 7-rays with energies greater 
than 20 MeV. The LAT has a peak effective area > 8000 cm^ at energies greater than 1 
GeV relevant for most of the event selections considered in this analysis and a large, ~ 2.4 
sr, field-of-view. The angular resolution for the 68% containment radius is ~ 0.6° for 1 
GeV photons that convert in the upper layers of the tracker [front events) and about a 
factor of 2 larger for those that convert in the bottom layers of the tracker (back events). A 
simple acceptance-averaged approximation for the 68% containment angle that is helpful to 
illustrate the energy dependent PSF is (^68(^)) = (0.8°) x (E/GeV) -"-^ ffi(0.07°). A full 
description of the LAT instrument is reported in lAtwood et al.l (120091 ). 



The data set used for the analysis of the AGNs includes LAT events with energy above 
100 MeV that were collected between 2008 August 4 and 2009 July 4. LAT-detected GRBs 
are considered up to 2009 September 30. A zenith angle cut of 105° was applied in order to 
greatly reduce any contamination from the Earth albedo. Blazars and GRBs have different 
emission characteristics, which result in different analysis procedures here. The event rate 
detected by the LAT in the vicinity (68% confidence radius) of a blazar is largely background 
dominated and only continuous observations over long time scales allow the detection of the 
underlying blazar emission. To minimize the background contamination when analyzing 
blazar data we use the "diffuse" class events, which provide the purest 7-ray sample and the 
best angular resolution. GRBs, on the other hand, emit most of their radiative 7-ray power 
on very short time scales (typically on the order of seconds) where the event rate can be 
considered mostly background free (at least during the prompt emission of bright bursts). 
It is therefore possible to loosen the event class selection to increase the effective area at 
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the expense of a higher background rate which is still small on short time scales for bright 
bursts. The "transient" class was designed for this specific purpose and we use these events 
for GRB analysis o 



2.1. AGN sample and potential biases 



We use blazars extracted from the First LAT AGN Catalog (ILAC; lAbdo et al.ll2010e[ ) 
as the AGN source sample to probe the UV through optical EBL. This catalog contains 671 
sources at high Galactic latitude (|6| > 10°) associated with high-confidence with blazars 
and other AGNs that were detected with a Test Statistiqj TS > 25 during the first 11 
months of science operation of the LAT. Detection of correlated multiwavelength variability 
was required in order to establish a source as being identified. 



Source associations were made with a Bayesian approach (similar to lMattox et al.ll200ll ). 
The Bayesian approach for source association implemented in the gtsrcid tool of the LAT Sci- 
enceToo/spackagqjuses only spatial coincidences between LAT and the counterpart sources . 
Candidate source catalogs u sed for this procedure in clude CRATES ([Healev et al.l 120071 ). 
CGRaBS faealev et aPboOsI ) and the Roma-BZCAT jMassaro et al.lboooh . which also pro- 
vide optical classifications and the latter two provide also spectroscopic redshifts for the 
sources. See lAbdo et al.l ( l2010ef ) for further details on the source detection and association 
algorithims refered to here. 

As discussed below, some methods applied here require one to distinguish among the 
different blazar source classes. Flat-spectrum radio quasars (FSRQs) and BL Lacs are dis- 
cerne d by their observed optical emission li ne equivalent widths and the Ca II break ratio 



(e.g..lStocke et al.lll99ll : lMarcha et al.lll996l ) following the procedure outlined in lAbdo et al. 



( l2009bl ). The BL Lac class itself is sub-divided into Low-, Intermediate-, and High- Synchrotron 
peaked BL Lacs (denoted as LSP-BLs, ISP-BLs and HSP-BLs, respectively) by estimating 
the position of their synchrotron peak, i^pg^k! from the indices aox — 0.384 • log(/5oooA//ikev) 
and ttro ^ 0.197- log(/5GHz//5000A) determined by the (rest frame) optical (/soooa), X-ray 
(/ikcv) a nd radio (/50Hz) flux densities listed in the online version of the Roma-BZCAT blazar 
catalog (jMassaro et al.ll2009l ). and using an empirical relationship between those broadband 



"see 



Atwood et al.l ( 20091 ) for further details on LAT event selection. 



'^The test statistic (TS) is defined as TS — —2 x {log{LQ) — log{Li)) with Lq the likelihood of the 
Null-hypothesis model as compared to the likelihood of a competitive model, Li, (see Section 3.2.2). 

^http://fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/overview. html 
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indices and i^p^ak ^^ derived in lAbdo et al.l ( l2010dl ). LSP-BLs have their synchrotron peak at 



K 



peak < 10^^ Hz, ISP-BLs at 10^^ Hz< z/^.^^ < 10^^ Hz and HS P-BLs at i^' y > 10 ^^ Hz. This 
is found to be in agreement with the classifications used in lAbdo et aLl (J2010d Jdl). Nearly 
all the 296 FSRQs are of LSP-type, while only 23% of the 300 BL Lacs are LSP-BLs, 15% 
are ISP-BLs, and 39% are HSP-BLs, 72 AGNs could not be classified, and 41 AGNs are of 
other type than listed above. 

Redshift information on the sources is extracted from the counterpart source catalogs 
(CRATES, CGRaBS, Roma-BZCAT). While all the redshifts of the ILAC FSRQs are known, 
only 42% of the high-confidence BL Lacs have measured redshifts. Obviously, AGN without 
redshift information are not used in the present work. 

The intrinsic average photon indices of Fermi blazars in the LAT energy range indicate 
a systematic hardening with source t ype from "^2.5 fo r FSRQs via ~ 2.2 for LSP- and 
~ 2.1 for ISP-, to ~ 1.9 for HSP-BLs JAbdo et al.lboiOch . On the other hand, their redshift 
distributions systematically decrease from the high- redshift (up to z ~ 3.1) FSRQs, via 
LSP-BLs located up to a redshift z ~ 1.5, down to the mostly nearby HSP-BLs at ^ < 0.5. 
This mimics a spectral softening with redshift if blazars are not treated separately by source 
type. A search for any systematic spectral evolution must therefore differentiate between 
the various AGN sub-classes (see below). 

To detect absorption features in the HE spectra (> 10 GeV; see Fig. [T]) of Fermi blazars, 
a thorough understanding of their intrinsic spectra, including variability, and source internal 
spectral features is required. Most blazars do not sh ow strong spectral variability in the 



LAT energy range o n > week scales (jAbdo et al.ll2010d ). despite often strong flux variability 



( lAbdo et al.ll2010gl ). Indeed at least three blazars which turn out to constrain the UV EBL 
the most (see Section 13.2.1(13.2.^ . show a > 99% probability of being variable in flux (using 
a x^ test) with a normalized excess variance of ~ 0.02 — 0.2 on time scales of hours to 
weeks. PKS 1502-1-106 (J1504-I-1029) is one of the most constraining sources in the sample. 
It displayed an exceptional flar e in August 2008 with a factor ~ 3 increase in flux within 
~ 12 hours (jAbdo et al.l 1201 Obi ). During this flare a flatter (when brighter) spectral shape 
was evident. The spectral curvature at the high energy end increased with decreasing flux 
level. If the high energy (> 10 GeV) photons are emitted during such flare activity, the 
constraints on the 7-ray optical depth would be tighter if only the flare-state spectral data 
were used. Because of limited photon statistics during the flare, however, we use the more 
conservative time averaged spectrum in the present analysis. 

Absorption in radiation fields internal to the source (e.g., accretion disk radiation, pho- 
ton emission from the broad line region) may cause a syste matic break in the 7-ray spectra 
that coincidentally mimics EBL attenuation (JReimerl 120071 ). In the case where such internal 
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absorption occurs, its redshift dependence is guaranteed, even in the absence of accretion 
evolution. T his is because of the redshifting of that energy where the interaction probabihty 



is maximum (iReimerl 120071 ) . Any technique that explores systematic variaton of observables 
(e.g., changes in spectral slope, flux ratios, e- folding cutoff energy) with redshift to single out 
EBL-induced absorption features in blazars with luminous accretion disk radiation (possibly 
indicated by strong emission lines) might therefore suffer from such a bias. 

All bright strong-line Fermi blazars (i.e., Fermz-FSRQs and some LSP- and ISP-BLs), 
howeve r, have been found to show spectral breaks already in the 1-10 GeV (source frame) 



range ( lAbdo et al.ll2010d ). This is too low in energy to be caused by EBL-attenuation for 
their redshift range < 3 (see Fig. [1]) . Although it is not clear if these breaks are due to internal 
absorption, the spectral softening results in low photon counts at energies > 10 GeV where 
EEL absorption is expected. Spectra of all bright HSP-BLs and some ISP-BLs, on the other 
hand, can be well represented by simple power laws without any signs of curvature or breaks. 
This indicates not only do they not have significant internal absorption in the 7-ray band, 
but also the absence of significant EBL absorption, which is expected to be beyond the LAT 
energy range for this nearby {z < 0.5) blazar population. 

Consequently, as we show in Section 13. H it remains challenging to quantify EBL ab- 
sorption effects in the LAT energy range based on population studies. On the other hand, 
the determination of the EBL-caused absorption features from individual blazars requires 
bright, high redshift objects with spectra extending to ^ lOGeV (Fig. [T]), and we focus on 
these blazars in Section 13.21 



2.2. GRB sample and potential biases 

The Fermi LAT has detected 11 GRBs from the beginning of its science operation 
(August, 4*^ 2008) until 30 September 2009, 6 of which have redshift measurements. Figure 
[3] shows the redshift and highest energy event associated with each of these GRBs. The 
probability of non-association is extremely small (see Table 1). 

GRB prompt emission is highly variable and shows signs of spectral evolution, a source 
of systematics to be considered carefully. Our approach in this paper is to restrict ourselves 
to the analysis of small time windows during the GRB emission where the temporal behavior 
does not seem to change significantly. 

The GRB spectral behavior is well-represented by the Band function (Band et al. 1993) 
in the keV-MeV range. An additional hard, F ~ 1.5-2, power-law comp onent, dominating 



at > 100 MeV, has now been firmly identified in a few GRBs: GRB090510 (jAckermann et al. 
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2010), GRB090902B, flAbdo et al.ll2009ch . GRB090926A flAbdo et al.ll2010hl ). Its absence 
in other LAT bursts could well be due to limited photon statistics. We assume that the 
power-law component extends well beyond 100 MeV up to ~ 10 GeV, below which EEL 
absorption is negligible (see Fig. [T]). EBL absorption is then expected to soften the power- 
law spectra from the extrapolation of the intrinsic/unabsorbed spectra beyond ~ 10 GeV. 
Systematic effects will, of course, occur when an intrinsic spectrum at high energies differs 
from this extrapolation. Source internal and/or intrinsic absorption via pair creation, e.g., 
would produce a curvature of the spectrum at higher energy which could be misinterpreted 
as an EBL absorption effect. Such a spectral break, whic h could be due to intrinsic pair 
creation, was detected in the LAT data from GRB 090926A (jAbdo et al.ll2010hl ) but we note 
that a corresponding roll-off in the intrinsic spectrum can only make our limits on the 7-ray 
optical depth more conservative. By contrast, a rising spectral component above >10 GeV 
would make our limits less constraining, but in the absence of any evidence for inverted 
gamma-ray spectra in GRBs, we consider this possibility unlikely. 



3. Analysis of 7-ray flux attenuation and results 

Assuming that high-energy photon absorption by the EBL is the sole mechanism that 
affects the 7-ray flux from a source at redshift 2;, the observed (i.e. absorbed) and unabsorbed 
fluxes at the observed energy E can be related by the optical depth, T^-y{E, z), as 

Fobs{E) = exp[-T^^{E, z)]Funabs{E). (1) 

This is the primary expression that we use to (i) explore 7-ray flux attenuation in the EBL 
from AGNs by means of a redshift-dependent flux ratio between a low- and a high- energy 
band; (ii) constrain EBL models which predict r^-y(-E, z) values much higher than the optical 
depth that would give the observed fluxes from individual blazars and GRBs; and (iii) put 
upper limits on the 7-ray optical depth calculated from the observed flux of individual blazars 
and GRBs, and the extrapolation of the unabsorbed flux to high energies. We discuss these 
methods and the results from our analysis below. 



3.1. Flux ratios - a population based method 



Because of inherent uncertainties in th e determinat i on of the intrinsic spectrum (Tint) 
for any given blazar in the pre- Fermi era, IChen et al.l (120041 ) proposed the average ratio 
F{> 10 GeV)/-F(> 1 GeV) for all blazars with signiflcant detections above 1 GeV, weighted 
according to the errors in -F(> 1 GeV), as a redshift-dependent tracer of the EBL attenuation 



-16- 



of 7-ray flux. The average flux ratio could then be compared with the predictions of the 
EBL models, taking selection effects into account. This approach assumes that the blazars 
are sampled from a homogeneous distribution with a single redshift-dependent luminosity 
function and a sing le intrinsic spectral index distribution. Preliminary results from Fermi 
( JAbdo et al.l l2009bl ) indicate that this assumption is inadequate. Consequently, we have 
calculated the average flux ratios for the different classes of blazars and discuss the results 
below. 

Among the AGN sample described in Section I^TT] we find that 237 FSRQs, 110 BL Lacs 
and 25 other AGNs are clearij ILAC associations with known redshift and detectable fluxes 
at energies > 1 GeV. There are 30 LSP-, 18 ISP- and 60 HSP- BL Lacs in this sub-sample. 

Of these AGN, only 22 FSRQs, 49 BL Lacs, and 1 other AGN have flux detections 
rather than upper limits above 10 GeV, including 10 LSP-, 6 ISP-, and 33 HSP-BL Lacs. 
For each of these BL Lacs and FSRQs, we calculated the ratio betw een the fluxes abo ve 
10 GeV and 1 GeV and their corresponding statistical errors following IChen et al.l (120041 ). 



Blazar type 


Num 


F 


ratio (pred) 


mean ratio (obs) 


red. x^ 


prob 


FSRQ 


22 


2.3 ±0.1 


0.04 ±0.01 


0.014 ±0.001 


4.38 


1.8 X 10-1° 


LSP-BL 


10 


2.2 ±0.1 


0.07 ±0.01 


0.028 ±0.003 


1.65 


0.11 


ISP-BL 


6 


2.1 ±0.1 


0.08 ±0.02 


0.048 ± 0.008 


1.86 


0.11 


HSP-BL 


33 


1.9 ±0.1 


0.12 ±0.03 


0.100 ±0.005 


1.29 


0.13 



Table 1: Spectral indices, mean predicted and observed flux ratios, and reduced x^ and 
probability for blazar sub-populations 

Figure |2] shows the observed flux ratios for the FSRQ population and BL Lac sub- 
populations as well as the ratios predicted according to the IFGL spectral index of each 
blazar, assuming an unbroken power law and no EBL attenuation. Table [T] shows the mean 
spectral index, mean flux ratios observed and expected, and the reduced x^ and associated 
probability given a parent distribution with constant flux ratio. As the blazar classes progress 
from FSRQ through LSP-BL, ISP-BL, and HSP-BL, 

1. the range of redshifts becomes narrower; 

2. on average, the spectra become harder; 



i.e., its association probability is at least 80%, it is the sole AGN associated with the correspond ing 



7-ray source, and it is not flagged to have problems that cast doubt on its detection ( Abdo et al.ll2010el ) 
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Fig. 2. — Flux ratio F{> 10 GeV)/F(> 1 GeV) as a function of redshift, in the Fermi LAT 
energy range, for FSRQs and BL Lac populations. The black diamonds are the observed 
ratios, while the red triangles show the ratio expected assuming an unbroken power law 
and no EBL attenuation. The black horizontal solid lines and red cross-hatched regions 
correspond to the mean observed ratios and expected ratios with errors. 
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3. both the predicted and observed mean flux ratios increase; 

4. the difference between the predicted and observed flux ratios decreases. 

The trend in the predicted flux ratios is a direct consequence of the hardening of the 
spectra as a function of source class, while the difference between the predicted and observed 
flux ratios is due to the fact that the curvature of the spectra decrease as the HE peak of 
the SED moves through the Fermi-LAT energies. The apparent discrepancies between the 
flux ratios for different blazar sub-populations arise from the fact that the LAT samples 
different parts of the blazar SED for these sub-classes. Indeed, a redshift distribution of the 
flux ratios for the combined blazar populations would show a strong, apparently decreasing 
trend, giving the appearance of an EBL absorption effect. When we separate the blazars 
into sub-populations, we find no significant redshift dependence of the flux ratios within 
each sub-population. The dearth of sources at high redshift and the large spread of spectral 
indices make it difficult to use the mean trend in the flux ratio as a function of redshift. 
To set upper limits on the 7-ray optical depth, we need to rely on the spectra of individual 
blazars, despite the increased dependence on the blazar emission model this entails. 

The flux ratio versus redshift relationship for BL Lacs is therefore primarily due to 
the differing intrinsic spectral characteristics of BL Lacs, rather than from EBL absorption. 
This test is a reminder of the importance of a careful consideration of the intrinsic spectral 
characteristics of the source populations chosen to probe EBL absorption. 



3.2. Constraints on EBL models from individual source spectra 

The sensitivity of the LAT over a broad energy range provides a unique opportunity to 
probe 7-ray spectra from AGNs and GRBs at < 10 GeV where EBL absorption is negligible 
and at > 10 GeV where EBL absorption can be substantial (see Fig. [T]). Thus extrapolations 
of the unabsorbed flux at low energies from individual sources to high energies, and assuming 
that the intrinsic spectra do not become harder at high energies, allows us to derive a measure 
of the total absorption (source in-situ and in EBL). We note that this is the only assumption 
made for the following methods. Furthermore, since any intrinsic spectral curvature or 
internal absorption effects can not be decoupled from EBL-caused curvature, the constraints 
derived below shall be considered as conservative upper limits on the EBL-caused opacity. 
These are then confronted with various EBL models. Clearly, high EBL density models 
possess a higher probability of being constrained by these methods than low density ones. 
In the following, we use two methods: the highest energy photon (Section 13.2. ip and the 
likelihood (Section I3.2.2p methods. 
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• BLLacs 

• FSRQs 
■ GRBs 



Kneiske - Best Fit 
Kneiske - High UV 
Salamon & Stecker - w corr. 
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Fig. 3. — Highest-energy photons from blazars and GRBs from different redshifts. Predic- 
tions of 77 optical depth r^^ = 1 (top panel) and r^^ = 3 (bottom panel) from various EBL 
models are indicated by lines. Photons above model predictions in this figure traverse an 
EBL medium with a high 7-ray opacity. The likelihood of detecting such photon considering 
the spectral characteristics of the source are considered in the method presented in section 

mm 
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Source 


Z 


Eraa-^ (GeV) 


Conv. Type 


AE/E 


68% radius 


Separation 


Chance Probability 


J1147-3812 


1.05 


73.7 


front 


10.7 % 


0.054° 


0.020° 


7.0 X 10""' 


(PKS 1144-379) 
















J1504+1029 


1.84 


48.9 


back 


5.4% 


0.114° 


0.087° 


5.6 X 10-3 


(PKS 1502+106) 




35.1 


back 


12.4% 


0.117° 


0.086° 


9.8 X 10^3 






23.2 


front 


7.2% 


0.072° 


0.052° 


5.6 X 10-3 


J0808-0751 


1.84 


46.8 


front 


9.7% 


0.057° 


0.020° 


1.5 X 10-3 


(PKS 0805-07) 




33.1 


front 


5.9% 


0.063° 


0.038° 


2.7 X 10-3 






20.6 


front 


8.9% 


0.075° 


0.029° 


6.9 X 10-3 


J1016+0513 


1.71 


43.3 


front 


11.4% 


0.054° 


0.017° 


1.2 X 10-3 


(CRATES J1016-I-0513) 




16.8 


front 


6.3% 


0.087° 


0.035° 


8.2 X 10-3 






16.1 


front 


7.6% 


0.084° 


0.018° 


8.2 X 10-3 


J0229-3643 


2.11 


31.9 


front 


10.7% 


0.060° 


0.035° 


1.7 X 10-3 


(PKS 0227-369) 
















GRB 090902B 


1.82 


33.4 


back 


10.5% 


0.117° 


0.077° 


6.0 X 10-* 


GRB 080916C 


4.24 


13.2 


back 


11.6% 


0.175° 


0.087° 


2.0 X 10-'' 



Table 2: List of blazars and GRBs detected by the LAT which have redshift measurements, 
and which constrain the EBL opacity the most. For each source, J2000 coordinate based 
name (other name), the energy of the highest-energy photon (HEP), the conversion type of 
the event {front or back) of the instrument, the energy resolution, AE/E, for 68% contain- 
ment of the reconstructed incoming photon energy, and the 68% containment radius based 
on the energy and incoming direction in instrument coordinates of the event, the separation 
from the source and the chance probability of the HEP being from the galactic diffuse or 
isotropic b ackgrounds are also liste d. The ener g y resol ution for the GRB HEP events is 
taken from lAbdo et al.l (l2009a ) and lAbdo et al.l (l2009dl ) using the respective lower energy 
bounds. The three highest-energy photons are listed for those sources that have multiple 
constraining photons. 
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3.2.1. Highest energy photons 

A simple method to constrain a given EBL model is to calculate the chance probability 
of detecting a photon with energy E > Emax, where Emax is the energy of the most energetic 
photon that we would expect when the source intrinsic spectrum is folded with the optical 
depth from the specific EBL model we want to test. We derive a conservative estimate of 
the intrinsic flux of the source by extrapolating the unabsorbed spectrum at low energies to 
high energies. We consider the LAT spectrum to be representative of the intrinsic spectrum 
at energies where the EBL is supposed to absorb less than ~ 1% of the photons for the most 
opaque models. This corresponds to an energy of around 10 GeV (down to ~ 6 GeV for GRB 
080916C at 2; ~ 4.3). Best fit spectral parameters of this "low-energy" unabsorbed spectrum 
were derived for all sources of the HEP set (see Table [3]). The spectrum is assumed to be 
a power law unless a significant deviation from this shape is measured at < 10 GeV (as is 
indeed observed from, e.g., FSRQs at GeV energies). This is the case for source J1504+1029, 
for which a log parabola model provides the best fit. 

Iterating through the source list described in Section 12.11 and Section 12.21 we find the 
energy -Emax of the highest-energy photon detected within the 68% containment radius (using 
the specific P6_V3_DIFFUSE instrument response functions for front and back events) of each 
source position. The resulting -Emax versus source redshift is shown in Figure [3] for sources 
with z > 0.5, and compared to the energy at which the optical depth r^^ is equal to 1 
and 3 according to the various EBL models. As shown in this figure, 5 AGN have -Emax 
that is sig nificantly grea t er (> 2) than the energy at which r^^ = 3 for the "baseline EBL 



model" of IStecker et al.l (120061 ). These 5 AGN (and 2 comparable GRBs) have emitted a 
number of events (hereafter highest- energy photons or HEP) that populate a region of the 
Emax — z phase space where EBL attenuation effects are predicted to be significant. These 
(henceforth called "HEP set" ) will be used in the following sections to constrain EBL models 
and to calculate the maximum amount of EBL attenuation that is consistent with the LAT 
observation^. 

It is possible that the high-energy photons considered here may not be emitted in the 
high-redshift source and instead are originated in any of the following background sources: 
Galactic 7-ray diffuse, isotropic (Extragalactic 7-ray plus charged-particle residuals) or a 
nearby point source. The likelihood of detecting any of these background photons within 
the 68% containment radius used to find the HEP set is quantified by computing the number 
of expected events within the 68% containment radius at the location of the source as deter- 



^Only the highest energy photon from each source is shown in Figured There are a few sources however 
with more than one constraining photon as indicated in Table [2j 
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mined by the best fit background model (Galactic and isotropic diffuse + point sources) and 
the instrument acceptance. The last column of Table [2] shows such probability for photons 
in the HEP set. These chance probabilities, although being fairly small, are non- negligible 
(at least in the case of blazars) if one would like to set significant constraints on specific 
EBL models by using this HEP. We later describe how this probability for the HEP to be 
a background fiuctuation was incorporated in our final results for this method. For now we 
will assume that this HEP is indeed from the source and we will first derive the type of 
constraints it allows us to set on different EBL models. We also note that a stricter set of 
cuts ( extrad iffuse) has been de veloped by the LAT team to study the Extragalactic 7-ray 



background (JAbdo et al.ll2010fl ). Despite the decreased 7-ray acceptance we find all photons 



in the HEP set to be retained when using these selection cuts. 

Monte-Carlo simulations are performed in order to test a particular EBL model with 
the derived intrinsic spectrum absorbed by the EBL as the Null-hypothesis. The simulations 
were performed using gtobssim, one of the science tools distributed by the Fermi science 
support center and the LAT instrument team. For each simulation we define the unabsorbed 
spectrum of the source as a power law (or log parabola in the case of J1504-I-1029) with 
spectral parameters drawn randomly from the best-fit values (and corresponding uncertainty) 
shown in Table [31 EBL absorption is applied according to the optical depth values of the 
considered model, and finally, the position and orientation of the Fermi satellite during the 
time interval described in Section 2 is used to account for the instrument acceptance that 
corresponds to the observations. The highest-energy photon from the simulated data is 
obtained following the exact same cuts and analysis procedure that was used for the data. 

The resulting distribution of the HEP simulated in each case (see e.g.. Figure S]) is 
used to estimate the chance probability of detecting a photon from the source with energy 
equal or greater than -Emax- We produced ~ 800000 and ~ 100000 simulations for each 
of the HEP sets for AGN and for GRBs, respectively. Assuming the HEP is indeed from 
the source, the probability of observing such high energy photon given the specific EBL 
model tested (called Phep) is calculated as the ratio between the number of cases where 
the HEP energy is above Emax (actually E^ax — <^Emax given the energy dispersion) and the 
total number of simulations performed. The number of simulations in each case was chosen 
to reach sufficient statistics at the tail of the distribution where the energy of the HEP is 
measured. Distributions of the HEP events from these MC simulations for GRB 080916C 
and GRB 090902B are shown in Figure [51 The open and filled histograms correspond to the 
distributi ons using the GRB s pectra without and with EBL absorption using the "baseline 



model" of IStecker et all tOOOJ ). 



To compute the final probability of rejection for the specific EBL model tested (called 
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Source 


normalization A'"o 
(10-^ph cm-2 s-i MeV-i) 


photon index F 
(PL) 


r, /3, E,/GeY 
(LP) 


TS 


J1147-3812 


0.570 ±0.081 


2.38 ±0.09 




221 


J1504+1029 


(1.84 ±0.23) -10-^ 




2.36 ±0.03, 

0.09 ±0.01, 

2.0 ±0.1 


34638 


J0808-0751 


1.212 ±0.078 


2.09 ±0.04 




1498 


J1016+0513 


1.183 ±0.078 


2.27 ±0.05 




1220 


J0229-3643 


0.789 ± 0.075 


2.56 ±0.07 




394 


J1012+2439 


0.552 ±0.058 


2.21 ±0.07 




443 


GRB 090902B 


146 ± 56 


1.40 ±0.37 




1956 


GRB 080916C 


1146 ±199 


2.15 ±0.22 




1398 



Table 3: Parameter values of the power law (PL) fits dN/dE = No{—T + 1)E ^ /[E^^^ — 
E-^+^] in the range E^^^ = 100 MeV to E^^^ = 10 GeV of the sources (AGN and GRBs) 
listed in Table H] except for source J1504+1029 where a log parabolic parametrization (LP) 
dN/dE = No{E / Eb)~^^~^f^^°^^^^^''^'^ has been found to be preferable over a power law fit 
(with ATS = 71.9). The spectral fits for the GRBs are performed below 6 GeV and 3 GeV 
for GRB 090902B and GRB 080916C, respectively. The TS values are obtained through a 
likelihood ratio test comparing a model with background only and a model where a point 
source was added. 
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Prejection), One needs to consider the fact that the HEP could be a background photon. We 
compute the probabihty for this to happen in Table 2 {Phkg)- In the end, one can fail to 
reject the EBL model because the HEP might be a background event or because there is a 
chance for a source photon with energy E^ax not to be absorbed by the EBL so that 

Prejection = Pbkg + PrEP X (1 ~ Pbkg) ■ (2) 

In Table 0, we list these 3 probabilities for each of our most constraining sources. When 
more than one photon is available for a given source, the probabilities are combined resulting 
in a stronger rejection. Although Prep can be quite constraining, our final significance of 
rejection is limited by Pbkg which is non-negligible in the case of blazars and which depends 
on the size of the region around each source defined a priori to look for associated high-energy 
events (68% PSF containment radius in this analysis). A larger HEP acceptance region (90% 
or 95% containment radius instead of 68%) would increase the background probability P^kg 
while also adding constraining photons to the HEP set. On a source- by-source basis, the 
rejection probability goes up or down with increasing radius depending on the number and 
energy of these additional photons, but our overall result remains the same. The unbinned 
likelihood method, which we describe in the next subsection (I3.2.2p . does not make use of 
an acceptance radius, and instead makes full use of available information in the data to 
systematically calculate a model rejection probability. 

The analysis described in t his section was app lied to all sources from the HEP set. 



We find the "baseline" model of IStecker et al.l ( l2006l ) to be significantly constrained by our 



observations. Colu mn 5 of Table H] shows the optical depth of the "baseli ne" model of 
Stecker et al.l (20061) for the HEP events. Since the "fast evolution" moderl of IStecker et al. 



( 120061 ) predicts higher opacities in the LAT energy range at all redshifts, our constraints on 



this model will naturally be higher than the ones found in Table H] for the "baseline" model. 



3.2.2. Likelihood method 

This second method to constrain specific EBL models makes use of a Likelihood-Ratio 
Test (LRT) technique. This approach compares the likelihood of the Null-hypothesis model 
(Lq) to best represent the data with the likelihood of a competitive model (-^i). The test 



^The "baseline" model considers the case where all galaxy 60/j,m luminosities evolved as (1 + z)^-^ up to 
z < 1.4, non-evolving between 1.4 < z < 6 and no emission at z > 6. In contrast, the "fast evolution" model 
assumes a more rapid galaxy luminosity evolution: ex (1 + z)"* for z < 0.8, oc (1 + z)^ for 0.8 < z < 1.5, no 
evolution for 1.5 < z < 6, and no emission at z > 6. Consequently, for a given redshift the "fast evolution" 
model predicts a higher 7-ray attentuation than the "baseline" model. 
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Source 


z 


-t^max 


TyZ, Emax) 


TyZ, -tLrnax) 


Number of photons 






(GeV) 


(F08) 


(St06, baseline) 


above 15 GeV 


Jl 147-3812 


1.05 


73.7 


0.40 


7.1 


1 


J1504+1029 


1.84 


48.9 


0.56 


12.2 


7 


J0808-0751 


1.84 


46.8 


0.52 


11.7 


6 


J1016+0513 


1.71 


43.3 


0.39 


9.0 


3 


J0229-3643 


2.11 


31.9 


0.38 


10.2 


1 


GRB 090902B 


1.82 


33.4 


0.28 


7.7 


1 


GRB 080916C 


4.24 


13.2 


0.08 


5.0 


1 



Table 4: Gamma-ray optical depth to HEP calculated usi np; the EEL ra o del o f Franceschini 
et al (2008; F08) in comparison to the "baseline" model of Istecker et~aD f l2006h (St06). Also 
listed are the number of photons associated to the source which have > 15 GeV energy and 
which can potentially constrain EBL models. 




lo Energy/GeV 



Fig. 4. — Distribution of highest-energy photons obtain ed from Mont e -Carlo simulations 
of the source J0808-0751 with the EBL attenuation by IStecker et al.l (|2006[ ). E^ax and 



E„ 



- (TEmax (wherc dEmax IS the energy uncertainty) are indicated by a solid and dotted 
vertical black lines, respectively. The probability of detecting a photon with energy equal or 
greater to Emax — (^Emax is equal to 6.8 x 10^^. 
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HEP method applied to Stecker 06 


HEP Rejection 


Source 


Z 


Energy (GoV) 


PbfeS 


Phep 


^rejection 


Significance 


J1147-3812 


1.05 


73.7 


7.0 X 10-* 


1.2 X 10-4 


8.1 X 10-* 


3.2 a 


J1504+1029 


1.84 


48.9 


5.6 X 10-3 


6.7 X 10-5 


5.7 X 10-3 








35.1 


9.8 X 10-3 


6.8 X 10-3 


1.7 X 10-2 








23.2 


5.6 X 10-3 


1.8 X 10-^ 


1.9 X 10-1 












Combined P,.ej = 1.7 X 10"^ 


4.1 o- 


J0808-0751 


1.84 


46.8 


1.5 X 10-3 


1.9 X 10-* 


1.7 X 10-3 








33.1 


2.7 X 10-3 


3.7 X 10-3 


6.4 X 10-3 








20.6 


6.9 X 10-3 


2.5 X 10-1 


2.6 X 10-1 












Combined P,.ej = 2.8 X lO-'' 


4.5 a 


J1016+0513 


1.71 


43.3 


1.1 X 10^3 


5.4 X 10-* 


1.6 X 10-3 








16.8 


8.2 X 10-3 


4.9 X IQ-i 


4.9 X 10-1 








16.1 


8.2 X 10-3 


6.5 X 10-1 


6.5 X 10-1 












Combined P^ej = 5.3 X IQ-* 


3.3o- 


J0229-3643 


2.11 


31.9 


1.7 X 10-3 


8.9 X 10-"^ 


1.8 X 10-3 


2.9 a 


GRB 090902B 


1.82 


33.4 


2 X lO-f" 


2.0 X 10-* 


2.0 X 10-* 


3.7 o- 


GRB 080916C 


4.24 


13.2 


8 X 10-** 


6.5 X 10-* 


6.5 X 10-* 


3.4 a 



Table 5: Listed are the significance of rejecting tlie "baseline" model (IStecker et al.l ( 120061 ) ). 
calculated using the HEP method as described in Section 13.2. 1[ For completeness, we also 
report individually the probability of the HEP to be a background event iPhkg) and the prob- 
ability for this HEP not to be absorbed by the EBL if it were emitted by the source (Phep) 
following Eq. |2l For those sources with more than one constraining photon, the individu al 
and combined PrpAprUnn. are calculated. The "fast evolution" model by lStecker et al.l (120061 ) is 



rejection 



more opaque and leads to an even higher significance of rejection. Applying this method to 
less opaque models leads to no hints of reject i on sin ce the probability Prep is large in those 
cases (e.g. > 0.1 for the iFranceschini et al.l (120081 ) EBL model). Note that a log parabola 
model was used as the intrinsic model for source J1504+1029 since evidence of curvature is 
observed here even below 10 GeV (see Table [3]). 
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statistic (TS) is defined as TS = — 2 x {log{Lo) — log{Li)). Following Wilks' theorem 
( IWilks Ill938[ ). the TS is asymptotically distributed as Xn (with n the difference in degrees of 



freedom betwee n the two models) if th e two models under consideration satisfy the following 
two conditions (jProtassov et al.l 120021 ): 1) the models must be nested and 2) the null- values 



of the additional parameters are not on the boundary of the set of possible parameter values. 

For the LRT we use the power law intrinsic spectrum convolved with the EBL absorption 
predicted by the model (Tmod) we are testing, exp[—aTmodiE,z)]Funabs{E), as the observed 
flux. For the Null-hypothesis we set a = 1 and we compare it to an alternative model where 
a is left as a free parameter, which therefore has one more degree of freedom than the Null- 
hypothesis. In the absense of any flux attenuation by the EBL, a = 0. Note that we allow 
the normalization parameter, a to go to negative values. This choice, although not physically 
motivated, allows us to satisfy the second condition mentioned above. As a consequence the 
Test-Statistics can simply be converted into a significance of rejecting the Null-hypothesis by 
making use of Wilks' theorem. Because of the lack of information on the intrinsic spectrum 
of a distant source above 10 GeV, we use the (unabsorbed) < 10 GeV observed spectrum as 
a reasonable assumption for the functional shape of the intrinsic source spectrum. A simple 
power law was found to be a good fit to the < 10 GeV data for the sources listed in Table H] 
except in the case of J1504+1029 where a log-parabolic spectrum was preferred. We note 
that if the actual intrinsic curvature is more pronounced than the one found with the best 
fit below 10 GeV, this would only make the results more constraining. 

As we mentioned earlier, although we are considering all EBL models in the literature, 
we find that our observations are only constraining the most opaque ones. Figure [6] shows 
the TS value as a function of the optical depth normalization parameter a, for the three 
most constraining blazars (J1016+0513, J0808-0751, J1504+1029) and the two GRBs ( GRBs 



090902B and 080916C) when considering the "baseline" model of IStecker et al.l (12006f ) with 
the LRT method. All sources are found to have an optical depth normalization parameter 
that is consistent with a > at the la level which is reassuring as we do not expect a 
rise in the spectrum on a physical basis. y/TSmax for a = 1 corresponds to the rejection 
significance for the specific model considered. The most constraining source, J1016+0513, 



rejects the Null- hypothesis (a = 1, corresponding to the "baseline" model of IStecker et al. 



( 120061 ) in this case) with a significance of ~ 6. Oct. This source could also constrain the "high 
UV model" of Kneiske et al. (2004) with a significance of 3.2(T although multi-trials effect 
substantially reduce this significance (see Section [3. 2. Sp . 

As compared to the HEP method, the LRT method incorporates the possibility of each 
photon being from the background into the unbinned maximum likelihood computation. 
Thus separate calculations of the background probablity and corresponding rejection proba- 
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10^ 




maximum energy (GeV) 



maximum energy (GeV) 



Fig. 5. — Distributions of the highest energy photons from simulations performed with esti- 
mates of our intrinsic spectra for GRB 080916C (l eft panel) and GRB 90902B (right panel), 
folded with EEL attenuation calculated using the IStecker et al.l (120061 ) baseline model. The 
total number of realizations (10^) in both the power-law and power-law convolved with the 
EBL cases is the same. 



Source 


z 


LRT Reje 
pre-trial 


'ction Significance 
post-trial 


J1147-3812 


1.05 


3.7cj 


2.0(7 


J1504+1029 


1.84 


4.6cj 


3.3 C7 


J0808-0751 


1.84 


5.4cj 


4.4 fT 


J1016+0513 


1.71 


e.Ocj 


5.1 a 


J0229-3643 


2.11 


3.2cj 


1.2 a 


GRB 090902B 


1.82 


3.6a 


1.9 a 


GRB 080916C 


4.24 


3.1cj 


1.0(7 



Table 6: Significance of rejecting the "baseline" model ( IStecker et al.l (120061 ) ). calculated 
using the LRT met hod described in Section 13.2.21 Again, the "fast evolution" model by 
Stecker et al.l (120061 ) leads to a high rejection significance with two sources (J0808-0751 and 
J1016+0513) with > 4o" post-trial significance. The post-trial significance is computed by 
taking into account the fact that our analysis is considering ~ 200 independent sources. 



-29- 



J1016+0513 (z=1.71) 
J0808-0751 (z=1.84) 
J1504+1029 (z=1.84) 
GRB 090902B {z=1.82) 
GRB 080916C (z=4.3) 




0.2 0.4 0.6 

Opacity normalization a 



Fig. 6. — Test statistic (TS) as a function of the 7-ray optical depth normahzation parameter 
calculated from the likelihood ratio test (LRT) for J1016 +0513, J0808- 751, J 1504+1029 and 
GRBs 090902B and 080916C. The "baseline" model of Istecker et"ZI tood ) has been used 
and the rejection for this model can be directly read out as ATS* between a = 1 and the best 
fit a for the source (horizontal dashed line). The confidence interval for the normalization 
parameter can be obtained using ATS* = CL^ where CL is the confidence level. 
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blity are not needed. Also since the LRT method takes into account all high-energy photons 
rather than the highest-energy ones in the HEP method, it gives more constraining results 
for the EBL model rejection with the exception of 2 GRBs where the HEP method gives 
slightly more constraining results. Finally, we note that the a priori choice of the size of the 
region around each source defined to look for associated high-energy events is a source of 
systematics for the HEP method ( which uses 68% PSF containment radius) while it does 
not affect the LRT method. 



3.2.3. Multi-trial effects and comhined probabilities 

Because the search for EBL signatures or rejection of specific EBL models is performed 
on all blazars and GRBs detected by the LAT, one has to consider multi-trials, which is 
potentially affecting our analysis. For independent searches, as is the case here, the post- 
trial probability threshold for obtaining a Aa result is -Ppost-triai = 1 — (1 — P4o-)^''^'"'''% where 
-^trials is the numbcr of trials and P4a- is the 4cr probability threshold for a s ingle search 
(^ 6.3 X 10^^). In the present case, the LAT AGN catalog that we have used fjAbdo et al. 
2010a) includes 709 AGNs of which ~ 200 have a sufficiently high redshift (~ 100 with 



> 10 GeV photon) to allow for the testing of EBL attenuation models with their 7-ray 
spectra. Only a handful of LAT GRBs were observed with sufficient statistics to hope to 
constrain the EBL. In the end, we have Nt^ais ~ 200 which corresponds to a post-trial 
probability for a Aa result of Pia,post-triai ~ 3.17 x 10~^. This corresponds to a significance 
of ~ 5. Her on an individual source which we will therefore consider as our threshold for a 
4cr post-trials rejection signi ficance for any specifi c EBL model. This Pia,post-triai threshold 
was reached in case of the IStecker et al.l (120061 ) "baseline" model for sources J0808-0751 
and J1016+0513 using the LRT method. Note that J1504+1029 is only slightly below this 
threshold. 

Combining specific EBL model rejection probabilities from multiple sourceqj we get a 
much hi gher rejection s i gnific ance. For HEP probabilities the combined rejection significance 
for the Istecker et all J2006h "baseline" model is ^ 8.9o- (^ 7.7a without the 2 GRBs) 
using Fishe r's method in order to combine results from independent tests of the same Null- 
hypothesis (JFisher Ill925l ). For the LRT method, we add the individual likelihood profiles 
to derive an overall profile from which y/TSmax gives an overall significance of 11.4cr for 
the same EBL model. Therefore both methods give very large reje c tion s ignificances even 
after taking multi-trial effects into account. Since the IStecker et al.l ( 120061 ) "fast-evolution" 



'since the spectral fits of all the sources we considered in this analysis are independent to each other. 
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model gives opacities larger than the "basehne model" in the LAT range, both models can 
be rejected by our analysis with very high confidence level. All other models can not be 
significantly rejected even after such stacking procedure is applied. 



3.3. Opacity upper limits 

Upper limits on the 7-ray optical depth have been evaluated with a method based on 
the comparison between the measured energy spectrum of the source and the unabsorbed 
spectrum above 10 GeV. The unabsorbed spectrum, Funabs, is assumed to be the extrapola- 
tion of the low-energy part, E<10 GeV, of the spectrum {Fe<io), where EBL attenuation is 
negligible (see Fig. [1]), to higher energies. F£;<io is fitted with a power-law or log-parabola 
function, according to the best TS value. At high energies, if no intrinsic hardening of the 
spectrum is present, the measured spectrum. Fobs, at (observed) energy E and the unab- 
sorbed spectrum, F^nabs, are related by Eq. [H The 7-ray optical depth can therefore be 
estimated at any given energy as 

T^^{E, z) = \n[F^r,abs{E)/Fo,,{E)]. (3) 

Since F^nabs is evaluated assuming no EBL attenuation, it gives a maximum value. Therefore 
the optical depth, t^^{E,z) given by Eq. [3] could already be considered as an upper limit, 
assuming that the difference between Funabs{E) and Fobs{E) is only due to EB L effects. The 



fit of both Fobs and F£;<io are carried out with a maximum likelihood analysis (JMattox et al. 



1996|) 



To evaluate -Fe;<io we have assumed a background model including all the point-like 
sources within 15° from the source under study and two diffuse components (Galactic and 
extra-galactic). The Galactic diffuse emission is modeled using a mapcube function, while a 
tabulated model is used for the isotropic component, representing the extragalactic emission 
as well as the residual instrumental background^ ^1. Both diffuse components are assigned a 
free normalization for the likelihood fit. In the fit we have considered all the nearby point 
sources within a 10° radius, modele d with a power-law with the photon index fixed to the 



value taken from the IFGL catalog (jAbdo et al.ll2010il ) and the integral fiux parameter left 



free. The remaining point sources are modeled with a power-law with all spectral parameters 
fixed. 



^^gtlike tool in the standard Fermi LAT Science Tools package provided by the Fermi Science Support 
Center (FSSC) 

^^http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
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The source under study has been fitted with a power-law and a log-parabola with all 
spectral parameters free. Among the two, we have chosen the fitted function showing the 
best TS value. The result is that for all the sources except J1504-I-1029 a power-law fit is 
preferred. From the fit results of -Fe<io we have extrapolated the spectral shape to obtain 
Funabs{E) above 10 GeV. 

A different method has been used to derive the measured fiux Fobs in selected energy 
bins. The whole energy range from 100 MeV to 100 GeV is divided in equal logarithmically 
spaced bins requiring in each energy bin a TS value greater than 10: 2 bins per decade above 
10 GeV for J0229-3643, J1016+0513 and J1147-3812, 4 bins per decade for J0808-0750 and 
5 bins per decade for J1504+1029. In each energy bin the standard gtlike tool has been 
applied assuming for all the point-like sources a simple power law spectrum with photon 
index fixed to 2.(o The integral fiuxes of all point-like sources within 10° are left as free 
parameters in the fitting procedure, while the diffuse background components are modeled 
as described in the previous paragraph. In this way, assuming that in each energy bin the 
spectral shape can be approximated with a power law, the flux of the source in all selected 
energy bins is evaluated. 

Once both Funabs and Fobs are determined, the maximum 7-ray optical depth in each 
energy bin can be estimated from Eq. (|3]). An upper limit on t^^{{E),z) with 95% CL 
in a constraining energy bin with mean energy (E) is then calculated by propagating the 
parameter uncertainties in the fltted fluJtj: 



r. 



-'y'y,UL95%CL{{E),z) = \n[Funabs{{E)) / Fobs{{E))] + 2a. (4) 



We compare these limits with the 7-ray optical depths predicted by various EBL models. 

In Figure [7] we show the upper limits (95% CL) derived at the mean energy of the bins 
above 10 GeV for various objects. In the highest energy bin the optical depth UL has been 
evaluated at the highest photon energy as reported in Table HI At this energy, the results 
of the optical depth UL at 99% CL are also reported (blue arrow). As an example, consider 
blazar J0808-0751 at z = 1.84 shown in the upper left plot: a larger optical depth would 
require an intrinsic spectrum that at high energies lies signiflcantly above the extrapolation 
obtained from the low energy spectrum. The flgure shows that the upper limit rules out 
those EBL models that predict strong attenuation. This result is consistent with all other 



^^since the energy bin is small enough to assume a flat spectrum. 

^■^It has been veryfied that the statistical errors follow a Gaussian distribution. The standard error prop- 
agation formula has then been applied. 
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upper limits derived with this method (see the other plots in Figure [7] and summary of the 
results in Table [7]). 



4. Discussion 

Studies with the highest energy extragalactic photons seen by the Fermi LAT primarily 
probe the UV and optical components of the EBL. The background fields responsible for 
7-ray attenuation can evolve strongly with redshift. In many of the models analyzed in this 
paper the EBL intensity can exceed the local value by a factor of 10 or more at redshifts 
near the peak of star-formation rate density. The optical depth to 7 rays from extragalactic 
sources is therefore determined by integrating the EBL intensity along the line of sight to 
the source from the observer. 

For an interaction angle of ^ = tt, the electron-positron pair production threshold con- 
dition leads to the value of the longest wavelength photon with which a 7-ray emitted at Zsrc 
with observed energy ii^obs can interact (source frame): 



Amax = 47.5 (1 + Zs 



E. 



obs 



GeV 



A. (5) 



The equation describes an upper limit on background photon wavelengths that can contribute 
to the 7-ray optical depth. Limits for Amax include 7175 A for blazar J1147-3812, and 4474 
A for GRB 090902B and 3286 A for GRB 080916C, based on the highest-energy 7-rays seen 
from these sources. In reality, interactions with shorter wavelength photons are more likely 
and will contribute more to the optical depth due to the redshifting of the 7-ray during 
propagation to Earth, the cross section, which peaks at approximately twice the threshold 
energy, and the geometry (interactions at angles of ~ 90° are more likely than head-on 
interactions as used in Eq. [5]). 

The results of our analysis of the highest energy 7-rays from blazars and GRBs detected 
by the Fermi LAT disfavor a UV ba ckground intensity a t the level predicted both by the 



baseline and fast-evolution models of IStecker et al.l (120061 ) . although the LAT observations 



discussed here do not constrain the predictions of this work at longer wavelengths. The 
two models of this work are based upon a backwards evolution model of galaxy formation. 
In this scenario, the IR SED of a galaxy is predicted from its luminosity at 60 fim. The 
locally-determined 60 /xm luminosity function is then assumed to undergo pure luminosity 
evolution following a power law in {1 + z). Optical and UV luminosities, relevant to Fermis 
extragalactic observations, are then determined by analytic approximation to the SEDs in 
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Fig. 7. — Derived upper limits for the optical depth of 7-rays emitted at z=1.84 (J0808-0751, 
J1504+1029), z=1.05 (J1147-3812) and z=1.71 (J1016+0513). Black arrows: upper limits 
at 95% c.l. in all energy bins used to determine the observed flux above 10 GeV. Red arrow: 
upper limits at 95% c.l. for the highest energy photon. Blue arrow: upper limit at 99% c.l. 
for the highest energy photon. The upper limits are inconsistent with the EBL models that 
predict the strongest opacity. 
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Salamon fc Steckerl (119981 ) and are normalized to the short wavelength portions of the IR 
SEDs. These models do not include absorption of UV light by dust in star-forming regions 
and the interstellar medium of galaxies, which may partially account for the high background 
in this model. While this model does account for redshift evolution in the UV-optical SEDs 
of galaxies, it does not allow for any evolut i on in the IR emission to which these SEDs are 
normalized. As mentioned by lStecker et al.l ( l2006l ). this is another factor which could result 
in overpredicted UV emission. 

Emissivity at UV wavelengths is closel y tied to the global star-formation rate density. 
Because the models of IStecker et al.l (120061 ) are not derived from an assumed function for 
the star-formation rate, limits on the UV emissivity in this case cannot be used to directly 
constrain star- formation. We do not find t hat o ur results are conclusively in disagreement 
with the 'best-fit' model of iKneiske et al.l (|2004[ ). In these models, the optical-UV EBL is 
based upon a Salpeter IMF and a star-formation rate density that peaks at z ~1.25, with 
a value of ~ 0.2 M© Mpc~^ yr~^) and falls slowly towards higher redshift. In the high-UV 
model, ultraviolet flux is boosted by a factor of 4 above the level of the best-fit model, greatly 
enhancing the opacity for 7-rays at energies below about 200 GeV. A star formation history 
of the magnitude required to produce the background in the high-UV m odel would be above 
essent ially all estimates of the global star formation rate (see for example iHopkins &: Beacom 
(I2OO6I )). All other EBL models are of such low density in the UV range that they can not 
be constrained by the data presented in this work. Although the results of our analysis can 
not yet place any stringent upper limits on the cosmological star-formation history that are 
competitive with current observational estimates, future prospects for probing low density 
UV models of the EBL by means of improved methods and enlarged GeV photon data sets 
may be promising. 

High-energy 7 rays that are absorbed by the EBL photons can initiate a pair cas- 
cade by subsequent Compton scattering of the CMB photons by the pairs. In case the 
intergalactic magnetic field (IGMF) is very weak, so that the pairs do not deflect out of 
our line of sight, this cascade radiation component can be detectable (|Plagalll995l). Ca l- 



culat i ons of such cascade signatures have been c arried out for AGNs (se e e.g. 



2OO2I: 



Murase et al 



2008 



Essev fc Kusenkd I2OIOI ) and for GRBs (see e.g. iDai &: Lu 



Dai et al. 



2002 



Razzaque et al.ll2004l : iTakahashi et al.ll2008l ) and found to compensate for a large portion of 



the flux that is abso rbed in the EBL. If blazars or GRBs are sources of ultrahigh energy cos- 
mic rays (UHECRs; IWaxman fc Coppil (119961 )). then photohadronic interactions by protons 
during their propagation in the background light can also induce a high-energy cascade sig- 
nature that would fo rm appreciable high-ene rgy emission, provided the IGM magnetic fleld 
is sufficiently small (JEssey fc Kusenkd l2010l ). However recent flux upper-limits calculated 
in the Fermi LAT range from TeV blazars lES 0347-121 and lES 0229+200 constrain the 
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IGMF to be > 3 X 10"^^ G (JNeronov fc Vovklboiol ). Such a strong field reduces the cascade 
flux significantly (because the emission becomes essentially isotropic due to the large deflec- 
tion angles) and the contribution to the observed flux is likely to be small. Furthermore, 
since the constraining blazar sample consists of FSRQs only, which seem weak TeV emitters, 
any of their reprocessed emission can only be small also. 

Exotic scenarios involving oscillation between 7-rays and axionlike particles, while prop- 
agating in the Galactic magneti c field, from distant sources may produce observable signa- 
tures in the TeV range (see e.g. iMirizzi et al.l 120071 : ISerpicd 120091 ). However the effect may 
not set in for typically assumed IGMF values or l ikely to be too small to ma ke up for EBL 
flux attenuation in the < 100 GeV range (see e.g. ISanchez-Conde et al.ll2009[ ) 



Conclusion 



Using the high-energy 11-month photon data set collected by Fermiirom distant blazars, 
and two GRBs we have (i) placed upper limits on the opacity of the Universe to 7 rays in the 
~10-100 GeV range coming from various redshifts up to 2; ~ 4.3; and ( ii) ruled out a i i EBL 
intensity in the redshift range ~ 1 to 4.3 as great as that predicted by lStecker et al.l (120061 ) 
in the ultraviolet range at more than 4a post trials in two independent sources (blazars). 
The overall rejection signiflcance is found to be > lOcr post trials therefore making this 
result very robust. Our most constraining sources are blazars J1504+1029, J0808-0751 and 
J1016+0513 with {z, (E^ax)) combinations of (1.84, 48.9 GeV), (1.84, 46.8 GeV) and (1.71, 
43.3 GeV), respectively. Although a likelihood ratio analysis of the latter source indicates 
that the sensitivity of our analysis method is approaching the EBL flux level of the "high UV 
model" of Kneiske et al (2004), multi-trial effects markedly reduced the rejection signiflcance. 
The two most constraining GRB s are GRB 090902B and GRB 080916C, both of which rule 
out the "baseline" EBL model of IStecker et al.l (12006 ) in the UV energy range at more than 
3cr level. The "fast evolution" model of IStecker et al.l (120061 ) predicts higher opacities in the 
LAT energy range at all reds hifts, and therefore is al s o ruled out. Together with the results 
from VHE ob s ervat ions (e.g., lAharonian et al.l (120071 ) : iMazin fc Raud (120071 )) the models by 
Stecker et al.l (120061 ) seem now disfavored in the UV and mid-IR energy range. We have also 
calculated model-independent optical depth upper-limits t^^,tjl{z, (Emax)) at 95% CL in the 
redshift 2; ~ 1 — 2.1 and -Emax ~ 28 — 74GeV ranges. 

As the high-energy photon data set collected by Fermi grows in the future and more 
blazars and GRBs are detected at constraining energies, the {E, z) phase space that con- 
strains r-y-y will become more populated. This will provide us with unique opportunities 
to constrain the opacity of the Universe to 7-rays over a large energy and redshift range. 
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and eventually help us further understand the evolution of the intensity of the extragalactic 
background light over cosmic time. 
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Source 


z 


-^max 


TUliz, Emax) 


Energy bins 
10 GeV-100 GeV 


J1147-3812 


1.05 


73.7 


1.33 


2 bins/dec 


J1504+1029 


1.84 


48.9 


1.82 


5 bins/dec 


J0808-0751 


1.84 


46.8 


2.03 


4 bins/dec 


J1016+0513 


1.71 


43.3 


0.83 


2 bins/dec 


J0229-3643 


2.11 


31.9 


0.97 


2 bins/dec 


J1012+2439 


1.81 


27.6 


2.41 


2 bins/dec 



Table 7: Upper limits (95% c.l.) on the 7-ray optical depth for AGN in Table ID The first 
and second column report the name of the sources and their redshift, the third column the 
maximum photon energy and the fourth column the optical depth UL evaluated at 95% c.l. 
as TuL = InlFunabsiE) / Fobs{E)] + 2a., the fourth column the number of energy bins/dec (for 
E>10 GeV) used to evaluate Fobs{E). 



